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Abstract 

A simple nondiagrammatic evaluation of the nondecoupling effect of heavy fermions 
on the Veltman's p parameter is presented in detail. This calculation is based on 
the path integral approach, the electroweak chiral Lagrangian formalism, and the 
Schwinger proper time method. 
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Although the Standard Model (SM) of particle physics was in very good agree- 
ment with the data during the past thirty years or so, there are strong indications 
that it is just a low-energy effective theory. For example, we do not discover the 
Higgs particle; we do not know how to unify the SM with gravity, why the SM gauge 
group is SU(3) C x SU(2)i x C/(l)y, and why there are 3 families of particles and so 
on; and the SM suffers the hierarchy problem, the unnaturalness problem, etc. At 
present there are many new physics extensions beyond the SM. Although we do not 
know whether nature really behaves like one of them or not, we can estimate their 
effects on the current electroweak precision observables. The Veltman's p parameter, 
defined by / (m^Cw) , is of particular interest among these observables p]. In the 
SM the p parameter is equal to 1 at tree level, which is protected by the custodial 
SU(2) symmetry in the Higgs sector of the SM. This symmetry overlaps with the 
local SU{2)l gauge symmetry. Mass splittings in weak isospin doublets break this 
symmetry. Their effect leaks into radiative corrections which make p differ slightly 
from 1. In model- independent studies the p parameter is related to the T parameter 
in the Peskin-Takeuchi's formalism [21 [3], and is also related to the /3i parameter in 
the electroweak chiral Lagrangian [U El El CH El El EE]- Whether the p parameter is 
1 or not is an indication whether the custodial symmetry is preserved or not in a 
generic theory. Current electroweak data fitting favors the p parameter being close 
to 1, that is, 1.00989 < p cxp < 1.01026 [EU [12], which puts stringent constraints 
on many new physics models. In literature the nondecoupling effects of a heavy 
fermion weak doublet or even an n-plet on the p parameter were already obtained 
using the Feynman diagrammatic calculation [31 HI [13]. In the present study we give 
a new calculational method of the contributions of heavy fermions to the p param- 
eter without referring to Feynman diagrams. The calculation itself is quite simple 
in the present method, although one needs some basic knowledge of the electroweak 
chiral Lagrangian and the Schwinger proper time method [TH [151 EEl ELY]- 

We start by considering a weak doublet of heavy fermions U and D, with masses 
■mu and m D respectively, and representations of SU{2) L x U(l)y as 

^=(p) Efc~(l,i), 2^(1,-1) (1) 

At some energy scale higher than the electroweak scale, the effective action of the 
heavy fermion sector can be written in a chiral invariant form as 
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-(ij L UMtlj R + ij R MU^ 



(2) 



where r a (a = 1,2,3) are Pauli matrices; g\ and #2 {Bp and W°) are the U{1)y and 
SU(2)l gauge couplings (fields), respectively; the dimensionless unitary unimodular 
matrix U (x) is the nonlinear realization of the Goldstone boson fields in the elec- 
troweak chiral Lagrangian; and the fermion mass matrix M = diag(mj/, m^). The 
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last two terms of eq. (J2]) will restore the usual fermion mass terms when the physical 
gauge, i.e. the unitary gauge U = 1, is taken. The U(x) field represents the de- 
gree of freedom of the would-be Goldstone bosons from the electroweak symmetry 
breaking, and its transformation under SU(2)l x U(1)y is given by 

U(x) -> V L (x)U(x)Vi(x) (3) 

where Vl(x) = exp{i^-9 a (x)} and Vr(x) = exp{i^-9°(x)} with 9 a (x) and 9°{x) being 
the SU{2) L and U(l)y group parameters, respectively. To derive the contribution 
of the heavy fermions to the p parameter, or equivalently the f3\ parameter in the 
electroweak chiral Lagrangian, we need to integrate out the heavy fermions above 
the electroweak scale which can be formulated as 

exp(z£W[^ W£, B,]) = J V^pVip exp(iS eS [U, W«, B„ i>}) (4) 

where Sew denotes the nondecoupling contribution of heavy fermions to the effective 
action just at the electroweak scale. Following Ref. [TJ, we can make the chiral 
decomposition of the U(x) field as follows: 

U(x)=Zl(x)Ux) (5) 

where £l{x) = exp{i : y0 a (x)} and £r(x) = exp{i^0°(a;)}, and their transformation 
under SU(2) L x U(1) Y are ' 

Cl(x) - h{x)£ L {x)Vl(x) , - h(x)Ux)V^x) (6) 

where h(x) = exp{i I ^-9h(x)} belongs to an induced hidden local U(l) symmetry 
group. Now we make a special SU(2)l x U(1)y chiral rotation as: Vl(x) = £l(x), 
Vr(x) = £r(x), under which U{x) is rotated to be 1 and the explicit ^-dependence 
of S e g disappears. Thus eq. (j2j) becomes 

Scs K, V* ^, ^} = j d 4 x^{i$ + f + ^ /5 -s- m)^ (7) 

where m = (mu+m£))/2, the traceless matrix s = Am-r 3 with Am = (m[/-mfl)/2, 
and the rotated fields are defined as 



^(x) = P L £ L {x)M*) + Pr£r(x)Mx) (8) 
1 r a r 3 

- 2 [92^v;{x) + gi - 



v,(x) = --[g/-V«(x) + 9l T -V«{x)] (9) 



a,(x) = -[ 92y y;(x) - g.-V^x)] (10) 

with ff2 ^;(x) = - ^]£i (n) 

3 3 

di^V^x) = ^[gJ-B^x) - id^ R (12) 
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According to eq.Q, the transformations of the above rotated fields under SU{2)i x 
U{1)y are given by 

i>t(x) - h(x)^(x) (13) 
92yV; - h(x)[g 2 ^ - id^x) (14) 

HijVS "> M*iy^° ~ i9 ^( x ) ( 15 ) 

Thus, the chiral symmetry SU(2)l x U(l)y covariance of the unrotated fields has 
been transferred totally to the hidden symmetry U(l) covariance of the rotated 
fields. We further find that a M (x) transforms covariantly: a M (x) — »■ h(x)a^(x)h'(x), 
while Vp(x) transforms as the "gauge field" of the hidden local U(l) symmetry: 
v ni x ) ~^ h{x)\v^{x)+idfj\}v{x). Accordingly, for any operator 0(x) which transforms 
covariantly under the hidden local symmetry: 0(x) — > h(x)0(x)hJ(x), we can define 
its covariant derivative as: d^O{x) = d fJi O{x) — i[v M (x), 0(x)]. It is important to 
note that the traces of some combinations of the rotated fields can be related to the 
terms of the electroweak chiral Lagrangian. First of all, from eqs. fllll) and (jl2j) it is 
straightforward to obtain: 

92^ K a (x) " g£v2(x) = -iUXA (16) 



where A M = U\D^U). Thus, the axial vector field a M = — ^rX^ r . And we can 
also write s in this covariant form: s = (fiAm • t 3 ^, since r 3 commutes with £r 
and £ R . On the other hand, we have 
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v, = ^rX^ r - i R9l T -B^ R + iUdA) (17) 



which leads to the following two relations: 



i.e. 



= -\xA-m^B,i R -{d^ R ) (is) 

1 r 3 
-iv^R = -CrX^ + z^^—B^- (d^ R ) (19) 



(d^ R ) + «W = (~\ X » - m^BJth (20) 
(d^ R ) - iv^ R = u\x, + igx^B^) (21) 



Thus, for any chiral rotated field / = £rF£ r , we have 



Ud»F)£ R + U\x, + i gi T -B,)F£} R + i R F{- l -X, - i 9l T -B^ R 
Z R {(D,F) + ~[X»,F\y R (22) 



where D^F = d^F + [igi^-B^, F]. In particular, if / = s = £i?Am ■ t 3 £ r , then 
F = Am ■ r 3 and we have 

d,s = ^|^(Am • r 3 ) + ±[X», Am. • r 3 ] |^ = ^Am^^, r 3 ]^, (23) 



and hence, 



d^s) = e^AmM,^] + ^Am[X„ [XV 3 ]] (24) 



Using the above relations, we can obtain the following identities: 

tr(sa M sa M ) = -^Am 2 tr(r 3 X M )tr(r 3 X M ) + ^Am 2 tr(X M X A 
tr[(cLs)(^ S )] = -tr[sdJd»s)) 



-Am 2 tr(r 3 X M )tr(r 3 X^) - Am 2 tr(X M A^) 



(25) 



(26) 



which are relevant to the /3i term in the electroweak chiral Lagrangian. Now, we 
proceed to integrate out the heavy fermion fields to get their contributions to the 
low-energy electroweak effective action. The integral measure of the heavy fermions 
remains unchanged under the special chiral rotation, since these fermions have an 
anomaly-free assignment of gauge charges. Thus, eq. (jl]) can be written as 



exp^Ewf^W;,^]) 



Vip^Dip^ exp 
VTp^Dip^ exp 



i / d A xip^(i$ + if) + ^75 — s — m)ip^ 



M 



d 4 xil>^($ — if> — z^7 5 — s + m)ijj£ 



(27) 



where the subscripts M and E indicate that the expressions are written respectively 
in the Minkowski spacetime and in the Euclidean spacetime (See Appendix [A] for 
the conversion relations of quantities in these two spaces). Here and henceforth, we 
mostly work in the Euclidean spacetime and, for convenience, we drop the subscript 
E in the expressions until further notice. Eq.( l271) gives 



iS EW [U,W^B ll }=\nDet(D 



m) 



Tr \\i{D + m) 



where D = $ — ii/i — ^75 — s. Since the imaginary part of the fermion determinant 
corresponds to the Wess-Zumino-Witten anomaly term [HI [191 [20] , f° r an Y anomaly- 
free underlying model we only need to consider the calculation of the real part of 
the fermion determinant. Using the Schwinger proper time formula [Hj, we have 



Re Tr \n(D + m) 



-Trln[(D t + m)(D + m)] 



— lim 

2 A^oo 



d A x I —tr cf i{x\e 
1 r 



r(D^+m)(D+m) 



X) (29) 



where tr c j^ denotes taking trace with respect to color (or technicolor), flavor, and 
Lorentz indices, and the operator (D' +m)(D + m) in the exponential can be further 
simplified as follows: 



+ m)(D + m) = E - V 2 + m 2 (30) 



with 



E = -2ms - 2zm^ 7s + - [7", 7"]^ + 7%« + ^{s, <vbs + s 2 (31) 

V M ee 9 M - - ia M 75 (32) 

Rpv = *[V M , = Vfa, + - 4a^)75 - i[a M > ( 33 ) 
ee - & - iv v ] = df,v v - d u Vf, - i[v M v u ] (34) 
ee d^s - i[v^ s] (35) 

Thus eq.f l29|) can be written as: 

ReTrlnp + m) = -- lim / rf 4 x / — e~ rm2 tr Cj/ii (a;|e- T(B " v2) |x) (36) 

2 A-^oo / / 1 T 

A 2 " 

where (x|e — ' r( -^ — v2 -* |x) is the so-called heat kernel, and it can be expanded in powers 
of r (and in powers of momenta as well) using the Seely-DeWitt expansion formula 
[15] as 



(x e 



-^Ix) = -^l\--+(-E 2 --[V II ,[W,E}}--R llu R^) 
1 1 16vr 2 \r 2 r V 2 6 L M L JJ 12 ^ ; 

+r( " + iV^' [V "' E]]E + E[V ^' [V "' E]] 



+ [V„E][^,E]))+^E^ + 0(r 3 



(37) 



From the structure of E and V M , we find that only the traces of the following terms 
in eq. fl37|) have contributions to the coefficient of tr (sa^sa^): 



tr, (^ E ^j ^ tea^ 



(3* 



T 



tr,( - ) 9 -16rm 2 sa M sa^ (39) 



tr/ ([V„ [V\ + £[V„ [V", E]} + [V M , £] [V", E}) \ 3 | , „. „„„ 



Tm 2 sa,,sa^ 



(40) 

tr,( — -E 4 ^) 9 ^T 2 m 4 S a^ (41) 
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where tr; stands for taking trace with respect to Lorentz spinor indices. Combing 
eqs ( IHSHTj) . we obtain 

tr^de-^-^lx) 3 -^(4 - — rm 2 + —^m^sa.sa^ (42) 
I67H 3 3 

which, by subsequently taking trace with respect to color (or technicolor) and flavor 
indices, leads to 

tr cJ j(x\e- TiE - V2) \x) 3 -^_(4 - ^rm 2 + ^-r 2 m A )ti f {sa^) (43) 

-LD7T o o 

where N c is the color (or technicolor) degree of freedom of the heavy fermions under 
consideration. From eq.( l26|) we see that tr[(<i M s)(<i M s)] and ti[sd^{d^s)] differ only 
by opposite sign, and thus we may recognize them as just one independent term. 
Likewise, by analyzing the structure of E and V M , we find that only the traces of 
the following terms in eq. (l3"T]) give contributions to the coefficient of tr[(rf M s)(<i Al s)]: 

trJ^E 2 ) 3 2(V)(d"a) (44) 

tai (j2 ( [V ^' [VM ' E]]E + E Fi» ^ + ^' E ^> 3 -^m 2 (d^s) 

(45) 

Combing the above equations gives 

tr z (x|e-^" v2 )|x) 3 JL(2 - ^rm 2 ){d,s){d^s) (46) 



which further implies 

tr cJ Ax\e-^ E -^\x) 3 ^(2 - ^rm 2 )tT f [(d,s)(d^s)] (47) 
Now, taking eqs. (1431) and (1471) into account, we finally arrive at: 

tr cJi/ (^|e" T(i? " V2) k) ^ T^((4-^rm 2 + ^r 2 m 4 )tr / (sa /1 ^) 

lD7T z [ O O 

+ (2-^rm 2 )tr / [(^)(# S )]} (48) 

In the following we will re-use the subscripts E and M to stand for, respectively, 
the quantities in the Euclidean space and those in the Minkowski space, and we will 
use tr as a short for try, i.e. taking trace over flavor space. Substituting eq. ()48i) into 
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eq. (l36|) . we obtain 

Re Tr hi(D + m) 

1 f * r dr _ Tm 2 N c l\ 40 2 16 2 4 . , u . 

3 -o / <Fx E hm / — e rm -— <^ (4 - — rm 2 + — r 2 m 4 tr (sa M sa M ) 

+ ( 2 -l rm 2 )tr[(^ S )(# S )] 
1 JV, 



2 16vr 2 



d 4 x s lim I [4(- 7 - In — ) - 8]tr(sa^) 

m 2 4 1 
[2(-7-ln^)-3]tr[(^)(d^)] 

d 4 x M A lim j [4(-7 -hi—)- 8]tr(sa IM sa tl ) 



2 16tt 2 



TO 2 4 1 

-[2(-7-ln^)-3]tr[(M(^)]j 



(49) 



where we have used eqs. ( 16211641) .( I56j) .(l57l) and f d A XE = if d A XM- Together with 
eqs.(|25D and (E51). eq.(H5j) gives 

Re Tr ln(D + m) 

i N, 



2 16tt 



1 rf 4 x M jim |[4(- 7 - ln^) - 8][-^Am 2 tr(r%)tr(r 3 ^)] 

?t? 2 4 1 1 
[2("7 - In^j) - -]-Am 2 tr(r%)tr(r 3 ^) j 

Am 2 y c/ 4 x Af tr(r 3 X At )tr(r 3 X' t ) (50) 



'24vr 2 



2 

where, as expected, the divergence term (—7 — In ^-J has been exactly canceled 
by the two parts. Thus, eq. fl50l) gives the nondecoupling contribution of the heavy 
fermions to the electroweak effective action. Comparing this term with the j3\ term 
of the standard electroweak chiral Lagrangian leads to, 

C'o = ^!/ 2 [tr(r 3 X,)] 2 = ^Am 2 [tr(r%)] 2 (51) 

which, together with the well-known relations, 4// 2 = e 2 /(s 2 c 2 M§) and Am = 
{mu — m£>)/2, implies, 

R - Nc e2 ( m u~ m p) 2 /„n 

Pl ~ 96vr 2 s 2 c 2 Ml 1 > 

where c = cos 9w, s = sin9w- Since p — 1 = a T = 2f3\, we finally obtain 

P ~ l + 487T 2 S 2 c 2 M\ [b6) 
8 



Eqs.f l52p and (1531) exactly coincide with the Feynman- diagram calculated results in 
previous works [3J HJ EES] • 

In summary, we have presented a new nondiagrammatic calculation method of 
the nondecoupling effect of heavy fermions on the p parameter. As we have seen, 
a crucial step of the present method is figuring out the coefficients of tr(sa M sa /i ) 
and tr[(d tl s)(d fl s)] in ReTrln(D + m). In the ordinary power counting of the QCD 
chiral Lagrangian, the scalar source term s is regarded as of order p 2 , and these two 
trace terms are of order p 6 and their coefficients were rarely considered in literature. 
However, in the power counting of the present case, s is just a constant, i.e. the 
mass splitting of the heavy fermions, and thus these two trace terms are of order 
p 2 and must give contributions to the (3\ term in the electroweak chiral Lagrangian. 
Our method is also applicable to the case of several heavy fermion doublets with 
mass mixing, and the generalization of this calculation would be straightforward. 



A Relations Between Minkowski Space and Eu- 
clidean Space 

To fix the notation, in this appendix we briefly review the relations between quan- 
tities in the Minkowski space and those in the Euclidean space. We use a mostly 
minus metric for the Minkowski space, {g^M — diag(l, —1, —I, —1), and we use a 
positive metric for the Euclidean space, (g^E = diag(l, 1, 1, 1). The coordinates in 
these two space are related by: 



-ix% , 



x 



M 



1,2,3) 



(54) 



which gives J d xm 



J d xe- From the correspondence p^ ~ d/(dx IJ '), we 



accordingly obtain the conversion relation of the momenta: (pq)m = i(Po)Ei (Pi)m = 
(Pi)e (i = 1, 2, 3), which by raising the indices gives, 

(p°)M = i(p°)E, (p 1 )m = -(p 1 )e (i = 1,2,3) (55) 

The vector source field v M (x) and the axial vector source field a M (x) are ordered to 
be converted exactly the same way as p M . Furthermore, a minus sign appears in 
the relation of the inner products of two vectors in these two spaces, {g^vP^P^M = 
(p°p° —p l p l )M = (—p°p° —p l p l )E = ~(9(iuP ,j 'P v )e, which suggests us to let the scalar 
source field s(x) be converted by: sm = —se- From these relations, we have 



tr(sa At sa A '^ 



-tr (sa^sa^ 



M 



tr[(d M s)(d M a)] 



-tr[(^s)(^s) 



(56) 
(57) 



On the other hand, since the Dirac matrices satisfy the Clifford algebra, {7^,7^} = 
2g^ u , with g^ LU different in these two spaces, their conversion relations are given by 



(7% = (W {l k )M = i{l k ) E (fe = l,2,3) 



(5* 
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In the Minkowski space, (7°)m is set to be hermitian, while the other 3 Dirac matrices 
are anti-hermitian. And thus (75 = "i7°7 1 7 2 7 3 )M is hermitian. As a result, in the 
Euclidean space all the Dirac matrices (j^e are hermitian, so is (75 = 7°7 1 7 2 7 3 )_b- 
Moreover, it is easy to check the relation: 

(7 5 )m = (7s)b (59) 
From eqs. (I5"5"|) and (I5"8|) . we have 



)M 



( 7 y - 7 y )m = i{iV + i k p% = m E (60) 



Likewise, ($)m = i($) E , (^)m = i{f)E, and (^)m = And we let the other 

quantities ip, ip, and m keep unchanged in these two spaces. Thus, finally we obtain 



cxp 
= exp 



i J d xip(i0 + f + ^75 — s — m) r (p 

d 4 xip($ — if — ijq§ — s + m)^ 



(61) 

J E 



B Necessary Integral Formulas 



In this appendix we list the integral formulas needed in the text as follows: 

lim / — e - rm = lim (-7 - In -^) (62) 

/OO 1 
(ire"™ 2 = — (63) 
1 m 2 

/OO 1 
drre~ Tm2 = — (64) 

where 7 is the Euler-Mascheroni constant, 7 ~ 0.5772. 
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